A regionally based precision medicine implementation initiative in North Africa:The PerMediNA consortium.

Precision Medicine is being increasingly used in the developed world to improve health care. While several Precision Medicine (PM) initiatives have been launched worldwide, their implementations have proven to be more challenging particularly in low- and middle-income countries. To address this issue, the "Personalized Medicine in North Africa" initiative (PerMediNA) was launched in three North African countries namely Tunisia, Algeria and Morocco. PerMediNA is coordinated by Institut Pasteur de Tunis together with the French Ministry for Europe and Foreign Affairs, with the support of Institut Pasteur in France. The project is carried out along with Institut Pasteur d'Algérie and Institut Pasteur du Maroc in collaboration with national and international leading institutions in the field of PM including Institut Gustave Roussy in Paris. PerMediNA aims to assess the readiness level of PM implementation in North Africa, to strengthen PM infrastructure, to provide workforce training, to generate genomic data on North African populations, to implement cost effective, affordable and sustainable genetic testing for cancer patients and to inform policy makers on how to translate research knowledge into health products and services. Gender equity and involvement of young scientists in this implementation process are other key goals of the PerMediNA project. In this paper, we are describing PerMediNA as the first PM implementation initiative in North Africa. Such initiatives contribute significantly in shortening existing health disparities and inequities between developed and developing countries and accelerate access to innovative treatments for global health.

A Longitudinal Study in Tunisia to Assess the Anti-RBD IgG and IgA Responses Induced by Three Different COVID-19 Vaccine Platforms.

BACKGROUND: Vaccination constitutes the best strategy against COVID-19. In Tunisia, seven vaccines standing for the three main platforms, namely RNA, viral vector, and inactivated vaccines, have been used to vaccinate the population at a large scale. This study aimed to assess, in our setting, the kinetics of vaccine-induced anti-RBD IgG and IgA antibody responses. METHODS: Using in-house developed and validated ELISA assays, we measured anti-RBD IgG and IgA serum antibodies in 186 vaccinated workers at the Institut Pasteur de Tunis over 12 months. RESULTS: We showed that RNA vaccines were the most immunogenic vaccines, as compared to alum-adjuvanted inactivated and viral-vector vaccines, either in SARS-CoV-2-naïve or in SARS-CoV-2-experienced individuals. In addition to the IgG antibodies, the vaccination elicited RBD-specific IgAs. Vaccinated individuals with prior SARS-CoV-2 infection exhibited more robust IgG and IgA antibody responses, as compared to SARS-CoV-2-naïve individuals. CONCLUSIONS: After following up for 12 months post-immunization, we concluded that the hierarchy between the platforms for anti-RBD antibody-titer dynamics was RNA vaccines, followed by viral-vector and alum-adjuvanted inactivated vaccines.

The African Human Microbiome Portal: a public web portal of curated metagenomic metadata.

There is growing evidence that comprehensive and harmonized metadata are fundamental for effective public data reusability. However, it is often challenging to extract accurate metadata from public repositories. Of particular concern is the metagenomic data related to African individuals, which often omit important information about the particular features of these populations. As part of a collaborative consortium, H3ABioNet, we created a web portal, namely the African Human Microbiome Portal (AHMP), exclusively dedicated to metadata related to African human microbiome samples. Metadata were collected from various public repositories prior to cleaning, curation and harmonization according to a pre-established guideline and using ontology terms. These metadata sets can be accessed at https://microbiome.h3abionet.org/. This web portal is open access and offers an interactive visualization of 14 889 records from 70 bioprojects associated with 72 peer reviewed research articles. It also offers the ability to download harmonized metadata according to the user's applied filters. The AHMP thereby supports metadata search and retrieve operations, facilitating, thus, access to relevant studies linked to the African Human microbiome. Database URL:  https://microbiome.h3abionet.org/.

Retrospective Phylodynamic and Phylogeographic Analysis of the Bluetongue Virus in Tunisia.

Bluetongue virus (BTV) is an arbovirus considered as a major threat for the global livestock economy. Since 1999, Tunisia has experienced several incursions of BTV, during which numerous cases of infection and mortality have been reported. However, the geographical origin and epidemiological characteristics of these incursions remained unclear. To understand the evolutionary history of BTV emergence in Tunisia, we extracted from Genbank the segment 6 sequences of 7 BTV strains isolated in Tunisia during the period 2000 to 2017 and blasted them to obtain a final dataset of 67 sequences. We subjected the dataset to a Bayesian phylogeography framework inferring geographical origin and serotype as phylodynamic models. Our results suggest that BTV-2 was first introduced in Tunisia in the 1960s and that since 1990s, the country has witnessed the emergence of other typical and atypical BTV serotypes notably BTV-1, BTV-3 and BTV-Y. The reported serotypes have a diverse geographical origin and have been transmitted to Tunisia from countries in the Mediterranean Basin. Interserotype reassortments have been identified among BTV-1, BTV-2 and BTV-Y. This study has provided new insights on the temporal and geographical origin of BTV in Tunisia, suggesting the contribution of animal trade and environment conditions in virus spread.

Longitudinal and Comparative Analysis of Gut Microbiota of Tunisian Newborns According to Delivery Mode.

Microbiota colonization is a dynamic process that impacts the health status during an individual's lifetime. The composition of the gut microbiota of newborns is conditioned by multiple factors, including the delivery mode (DM). Nonetheless, the DM's influence remains uncertain and is still the subject of debate. In this context, the medical indication and the emergency of a cesarean delivery might have led to confounding conclusions regarding the composition and diversity of the neonatal microbiome. Herein, we used high-resolution shotgun sequencing to decipher the composition and dynamics of the gut microbiota composition of Tunisian newborns. Stool samples were collected from 5 elective cesarean section (ECS) and 5 vaginally delivered (VD) newborns at the following time points: Day 0, Day 15, and Day 30. The ECS and VD newborns showed the same level of bacterial richness and diversity. In addition, our data pointed to a shift in microbiota community composition during the first 2 weeks, regardless of the DM. Both ECS and VD showed a profile dominated by Proteobacteria, Actinobacteria, and Firmicutes. However, ECS showed an underrepresentation of Bacteroides and an enrichment of opportunistic pathogenic species of the ESKAPE group, starting from the second week. Besides revealing the intestinal microbiota of Tunisian newborns, this study provides novel insights into the microbiota perturbations caused by ECS.

Phylogenetic and amino acid signature analysis of the SARS-CoV-2s lineages circulating in Tunisia.

Since the beginning of the Coronavirus disease-2019 pandemic, there has been a growing interest in exploring SARS-CoV-2 genetic variation to understand the origin and spread of the pandemic, improve diagnostic methods and develop the appropriate vaccines. The objective of this study was to identify the SARS-CoV-2s lineages circulating in Tunisia and to explore their amino acid signature in order to follow their genome dynamics. Whole genome sequencing and genetic analyses of fifty-eight SARS-CoV-2 samples collected during one-year between March 2020 and March 2021 from the National Influenza Center were performed using three sampling strategies.. Multiple lineage introductions were noted during the initial phase of the pandemic, including B.4, B.1.1, B.1.428.2, B.1.540 and B.1.1.189. Subsequently, lineages B1.160 (24.2%) and B1.177 (22.4%) were dominant throughout the year. The Alpha variant (B.1.1.7 lineage) was identified in February 2021 and firstly observed in the center of our country. In addition, A clear diversity of lineages was observed in the North of the country. A total of 335 mutations including 10 deletions were found. The SARS-CoV-2 proteins ORF1ab, Spike, ORF3a, and Nucleocapsid were observed as mutation hotspots with a mutation frequency exceeding 20%. The 2 most frequent mutations, D614G in S protein and P314L in Nsp12 appeared simultaneously and are often associated with increased viral infectivity. Interestingly, deletions in coding regions causing consequent deletions of amino acids and frame shifts were identified in NSP3, NSP6, S, E, ORF7a, ORF8 and N proteins. These findings contribute to define the COVID-19 outbreak in Tunisia. Despite the country's limited resources, surveillance of SARS-CoV-2 genomic variation should be continued to control the occurrence of new variants.

The Delta variant wave in Tunisia: Genetic diversity, spatio-temporal distribution and evidence of the spread of a divergent AY.122 sub-lineage.

Introduction: The Delta variant posed an increased risk to global public health and rapidly replaced the pre-existent variants worldwide. In this study, the genetic diversity and the spatio-temporal dynamics of 662 SARS-CoV2 genomes obtained during the Delta wave across Tunisia were investigated. Methods: Viral whole genome and partial S-segment sequencing was performed using Illumina and Sanger platforms, respectively and lineage assignemnt was assessed using Pangolin version 1.2.4 and scorpio version 3.4.X. Phylogenetic and phylogeographic analyses were achieved using IQ-Tree and Beast programs. Results: The age distribution of the infected cases showed a large peak between 25 to 50 years. Twelve Delta sub-lineages were detected nation-wide with AY.122 being the predominant variant representing 94.6% of sequences. AY.122 sequences were highly related and shared the amino-acid change ORF1a:A498V, the synonymous mutations 2746T>C, 3037C>T, 8986C>T, 11332A>G in ORF1a and 23683C>T in the S gene with respect to the Wuhan reference genome (NC_045512.2). Spatio-temporal analysis indicates that the larger cities of Nabeul, Tunis and Kairouan constituted epicenters for the AY.122 sub-lineage and subsequent dispersion to the rest of the country. Discussion: This study adds more knowledge about the Delta variant and sub-variants distribution worldwide by documenting genomic and epidemiological data from Tunisia, a North African region. Such results may be helpful to the understanding of future COVID-19 waves and variants.

First Detection of Human ST131-CTX-M-15-O25-B2 Clone and High-Risk Clonal Lineages of ESBL/pAmpC-Producing E. coli Isolates from Diarrheic Poultry in Tunisia.

Circulation of a multi-resistance clone of bacteria associated with genetic elements in diseased animals constitutes a global public health problem. Our study focused on the characterization of the support of ESBL in cefotaxime resistant E. coli (CTXR) isolates recovered from poultry with diarrhea, analysis of their clonal lineage, and virulence-associated genes. The study was carried out on 130 samples of chickens with diarrhea, collected in 2015 from poultry farms in Tunisia. Isolates of 20 CTXR E. coli strains were identified as ESBL and AmpC ß- lactamase producers. The following ß-lactamase genes (number of isolates) were detected: blaCTX-M-15+ blaOXA1 (4), blaCTX-M-15 + blaOXA1 + blaTEM-1b (2), blaCTX-M-1 + blaTEM-1b (9), blaCTX-M-1 (2), blaCMY2 + blaTEM-1b (3). Six E. coli harboring blaCTXM-15 were allocated to ST131-B2-O25b-; six and three blaCTX-M-1 were grouped in ST155, ST10, and ST58, respectively, related to the phylogroup D and A. The qnrB gene, the variant aac(6')-Ib-cr, and the class 1 integrons with different gene cassettes, were detected amongst our 20 isolated strains, which were classified as ExPEC and aEPEC. Our findings highlighted the emergence of the human pandemic ST131-CTX-M-15-O25-B2 clone and the high risk of such clonal lineage strains in diarrheic poultry, in Tunisia, which could constitute a risk of their transfer to healthy animals and humans.

Human OMICs and Computational Biology Research in Africa: Current Challenges and Prospects.

Following the publication of the first human genome, OMICs research, including genomics, transcriptomics, proteomics, and metagenomics, has been on the rise. OMICs studies revealed the complex genetic diversity among human populations and challenged our understandings of genotype-phenotype correlations. Africa, being the cradle of the first modern humans, is distinguished by a large genetic diversity within its populations and rich ethnolinguistic history. However, the available human OMICs tools and databases are not representative of this diversity, therefore creating significant gaps in biomedical research. African scientists, students, and publics are among the key contributors to OMICs systems science. This expert review examines the pressing issues in human OMICs research, education, and development in Africa, as seen through a lens of computational biology, public health relevant technology innovation, critically-informed science governance, and how best to harness OMICs data to benefit health and societies in Africa and beyond. We underscore the disparities between North and Sub-Saharan Africa at different levels. A harmonized African ethnolinguistic classification would help address annotation challenges associated with population diversity. Finally, building on the existing strategic research initiatives, such as the H3Africa and H3ABioNet Consortia, we highly recommend addressing large-scale multidisciplinary research challenges, strengthening research collaborations and knowledge transfer, and enhancing the ability of African researchers to influence and shape national and international research, policy, and funding agendas. This article and analysis contribute to a deeper understanding of past and current challenges in the African OMICs innovation ecosystem, while also offering foresight on future innovation trajectories.

H3ABioNet genomic medicine and microbiome data portals hackathon proceedings.

African genomic medicine and microbiome datasets are usually not well characterized in terms of their origin, making it difficult to find and extract data for specific African ethnic groups or even countries. The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for developing data portals for African genomic medicine and African microbiomes to address this and ran a hackathon to initiate their development. The two portals were designed and significant progress was made in their development during the hackathon. All the participants worked in a very synergistic and collaborative atmosphere in order to achieve the hackathon's goals. The participants were divided into content and technical teams and worked over a period of 6 days. In response to one of the survey questions of what the participants liked the most during the hackathon, 55% of the hackathon participants highlighted the familial and friendly atmosphere, the team work and the diversity of team members and their expertise. This paper describes the preparations for the portals hackathon and the interaction between the participants and reflects upon the lessons learned about its impact on successfully developing the two data portals as well as building scientific expertise of younger African researchers. Database URL: The code for developing the two portals was made publicly available in GitHub repositories: [https://github.com/codemeleon/Database; https://github.com/codemeleon/AfricanMicrobiomePortal].

Phylogenetic analysis of the neuraminidase segment gene of Influenza A/H1N1 strains isolated from Monastir Region (Tunisia) during the 2017-2018 outbreak.

Influenza A/H1N1 is widely considered to be a very evolutionary virus causing major public health problems. Since the pandemic of 2009, there has been a rapid rise in human Influenza virus characterization. However, little data is available in Tunisia regarding its genetic evolution. In light of this fact, our paper aim is to genetically characterize the Neuraminidase, known as the target of antiviral inhibitors, in Tunisian isolates circulating in Monastir region during 2017-2018. In total of 31 positive Influenza A/H1N1 detected by multiplex real-time PCR, RT-PCR of neuraminidase was performed. Among the 31 positive samples, 7 samples representing fatal and most severe cases were conducted for sequencing and genetic analysis. The results thus obtained showed genetic evolution of the A/H1N1 neuraminidase between 2009 and 2010 and 2018-2019 outbreaks. All Tunisian isolates were genetically related to the recommended vaccine strain with a specific evolution. Moreover, the phylogenetic analysis demonstrated that France and especially Italian strains were the major related strains. Interestingly, our results revealed a specific cluster of Tunisian isolates where two intragroup were evolved in correlation with the severity and the fatalities cases. From the outcome of our investigation, this study confirms the genetic evolution of the Influenza A virus circulating in Tunisia and gives a preliminary analysis for a better comprehension of new emerging Tunisian strain's virulence and thus, a more appropriate monitoring of Influenza virus A/H1N1 during each round of outbreaks. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11756-021-00723-y.

Historical origins and zoonotic potential of avian influenza virus H9N2 in Tunisia revealed by Bayesian analysis and molecular characterization.

During 2009-2012, several outbreaks of avian influenza virus H9N2 were reported in Tunisian poultry. The circulating strains carried in their hemagglutinins the human-like marker 226L, which is known to be important for avian-to-human viral transmission. To investigate the origins and zoonotic potential of the Tunisian H9N2 viruses, five new isolates were identified during 2012-2016 and their whole genomes were sequenced. Bayesian-based phylogeny showed that the HA, NA, M and NP segments belong to the G1-like lineage. The PB1, PB2, PA and NS segments appeared to have undergone multiple intersubtype reassortments and to be only distantly related to all of the Eurasian lineages (G1-like, Y280-like and Korean-like). The spatiotemporal dynamic of virus spread revealed that the H9N2 virus was transferred to Tunisia from the UAE through Asian and European pathways. As indicated by Bayesian analysis of host traits, ducks and terrestrial birds played an important role in virus transmission to Tunisia. The subtype phylodynamics showed that the history of the PB1 and PB2 segments was marked by intersubtype reassortments with H4N6, H10N4 and H2N2 subtypes. Most of these transitions between locations, hosts and subtypes were statistically supported (BF > 3) and not influenced by sampling bias. Evidence of genetic evolution was observed in the predicted amino acid sequences of the viral proteins of recent Tunisian H9N2 viruses, which were characterized by the acquisition of new mutations involved in virus adaptation to avian and mammalian hosts and amantadine resistance. This study is the first comprehensive analysis of the evolutionary history of Tunisian H9N2 viruses and highlights the zoonotic risk associated with their circulation in poultry, indicating the need for continuous surveillance of their molecular evolution.

Investigation of serum proteome homeostasis during radiation therapy by a quantitative proteomics approach.

The purpose of the present study is to analyze the serum proteome of patients receiving Radiation Therapy (RT) at different stages of their treatment to discovery candidate biomarkers of the radiation-induced skin lesions and the molecular pathways underlying the radiation signatures. Six stages of RT treatment were monitored from patients treated because of brain cancer: before starting the treatment, during the treatment (four time points), and at 4 weeks from the last RT dose. Serum samples were analyzed by a proteomics approach based on the Data Independent Acquisition (DIA) mass spectrometry (MS). RT induced clear changes in the expression levels of 36 serum proteins. Among these, 25 proteins were down- or up-regulated significantly before the emergence of skin lesions. Some of these were still deregulated after the completion of the treatment. Few days before the appearance of the skin lesions, the levels of some proteins involved in the wound healing processes were down-regulated. The pathway analysis indicated that after partial body irradiation, the expression levels of proteins functionally involved in the acute inflammatory and immune response, lipoprotein process and blood coagulation, were deregulated.

Association of common IL-10 promoter gene variants with the susceptibility to head and neck cancer in Tunisia

Background/aim: We investigated the association of three IL-10 promoter single-nucleotide polymorphisms and altered IL-10 plasma levels with the risk of head and neck cancer (HNC). Materials and methods: Study subjects comprised 194 HNC patients [137 nasopharyngeal cancer (NPC) and 57 laryngeal cancer (LC)], and 263 healthy controls. Genotyping of rs1800896 (-1082A>G), rs1800871 (-819C>T), and rs1800872 (-592A>C) IL-10 variants was performed by real-time PCR; IL-10 levels were measured by enzyme amplified immuno sensitivity assay (EAISA). Results: Study subjects comprised 194 HNC patients [137 nasopharyngeal cancer (NPC) and 57 laryngeal cancer (LC)], and 263 healthy controls. Genotyping of rs1800896 (-1082A>G), rs1800871 (-819C>T), and rs1800872 (-592A>C) IL-10 variants was performed by real-time PCR; IL-10 levels were measured by enzyme amplified immuno sensitivity assay (EAISA). Conclusion: Our results demonstrate that IL-10-1082, IL-10-819, and IL-10-592 variants, and haplotypes GC and GT constitute biomarkers for early detection of HNC, especially NPC subtype. IL-10 -819T/C and TA haplotype may be used as biomarkers for early detection of LC.

Developing reproducible bioinformatics analysis workflows for heterogeneous computing environments to support African genomics.

BACKGROUND: The Pan-African bioinformatics network, H3ABioNet, comprises 27 research institutions in 17 African countries. H3ABioNet is part of the Human Health and Heredity in Africa program (H3Africa), an African-led research consortium funded by the US National Institutes of Health and the UK Wellcome Trust, aimed at using genomics to study and improve the health of Africans. A key role of H3ABioNet is to support H3Africa projects by building bioinformatics infrastructure such as portable and reproducible bioinformatics workflows for use on heterogeneous African computing environments. Processing and analysis of genomic data is an example of a big data application requiring complex interdependent data analysis workflows. Such bioinformatics workflows take the primary and secondary input data through several computationally-intensive processing steps using different software packages, where some of the outputs form inputs for other steps. Implementing scalable, reproducible, portable and easy-to-use workflows is particularly challenging. RESULTS: H3ABioNet has built four workflows to support (1) the calling of variants from high-throughput sequencing data; (2) the analysis of microbial populations from 16S rDNA sequence data; (3) genotyping and genome-wide association studies; and (4) single nucleotide polymorphism imputation. A week-long hackathon was organized in August 2016 with participants from six African bioinformatics groups, and US and European collaborators. Two of the workflows are built using the Common Workflow Language framework (CWL) and two using Nextflow. All the workflows are containerized for improved portability and reproducibility using Docker, and are publicly available for use by members of the H3Africa consortium and the international research community. CONCLUSION: The H3ABioNet workflows have been implemented in view of offering ease of use for the end user and high levels of reproducibility and portability, all while following modern state of the art bioinformatics data processing protocols. The H3ABioNet workflows will service the H3Africa consortium projects and are currently in use. All four workflows are also publicly available for research scientists worldwide to use and adapt for their respective needs. The H3ABioNet workflows will help develop bioinformatics capacity and assist genomics research within Africa and serve to increase the scientific output of H3Africa and its Pan-African Bioinformatics Network.

Organizing and running bioinformatics hackathons within Africa: The H3ABioNet cloud computing experience.

The need for portable and reproducible genomics analysis pipelines is growing globally as well as in Africa, especially with the growth of collaborative projects like the Human Health and Heredity in Africa Consortium (H3Africa). The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for portable, reproducible pipelines adapted to heterogeneous compute environments, and for the nurturing of technical expertise in workflow languages and containerization technologies. To address this need, in 2016 H3ABioNet arranged its first Cloud Computing and Reproducible Workflows Hackathon, with the purpose of building key genomics analysis pipelines able to run on heterogeneous computing environments and meeting the needs of H3Africa research projects. This paper describes the preparations for this hackathon and reflects upon the lessons learned about its impact on building the technical and scientific expertise of African researchers. The workflows developed were made publicly available in GitHub repositories and deposited as container images on quay.io.

Phagonaute: A web-based interface for phage synteny browsing and protein function prediction.

Distant homology search tools are of great help to predict viral protein functions. However, due to the lack of profile databases dedicated to viruses, they can lack sensitivity. We constructed HMM profiles for more than 80,000 proteins from both phages and archaeal viruses, and performed all pairwise comparisons with HHsearch program. The whole resulting database can be explored through a user-friendly "Phagonaute" interface to help predict functions. Results are displayed together with their genetic context, to strengthen inferences based on remote homology. Beyond function prediction, this tool permits detections of co-occurrences, often indicative of proteins completing a task together, and observation of conserved patterns across large evolutionary distances. As a test, Herpes simplex virus I was added to Phagonaute, and 25% of its proteome matched to bacterial or archaeal viral protein counterparts. Phagonaute should therefore help virologists in their quest for protein functions and evolutionary relationships.

Detection of generic differential RNA processing events from RNA-seq data.

RNA-seq data analysis has revealed abundant alternative splicing in eukaryotic mRNAs. However, splicing is only one of many processing events that transcripts may undergo during their lifetime. We present here RNAprof (RNA profile analysis), a program for the detection of differential processing events from the comparison of RNA-seq experiments. RNAprof implements a specific gene-level normalization procedure and compares RNA-seq coverage profiles at nucleotide resolution to detect regions of significant coverage differences, independently of splice sites or other gene features. We used RNAprof to analyze the effect of alternative-splicing regulators NSRa and NSRb on the Arabidopsis thaliana transcriptome. A number of intron retention events and alternative transcript structures were specifically detected by RNAprof and confirmed by qRT-PCR. Further tests using a public Mus musculus RNA-seq dataset and comparisons with other RNA isoform predictors showed that RNAprof uniquely identified sets of highly significant processing events as well as other relevant library-specific differences in RNA-seq profiles. This highlights an important layer of variation that remains undetected by current protocols for RNA-seq analysis.

In silico prediction of protein-protein interactions in human macrophages.

BACKGROUND: Protein-protein interaction (PPI) network analyses are highly valuable in deciphering and understanding the intricate organisation of cellular functions. Nevertheless, the majority of available protein-protein interaction networks are context-less, i.e. without any reference to the spatial, temporal or physiological conditions in which the interactions may occur. In this work, we are proposing a protocol to infer the most likely protein-protein interaction (PPI) network in human macrophages. RESULTS: We integrated the PPI dataset from the Agile Protein Interaction DataAnalyzer (APID) with different meta-data to infer a contextualized macrophage-specific interactome using a combination of statistical methods. The obtained interactome is enriched in experimentally verified interactions and in proteins involved in macrophage-related biological processes (i.e. immune response activation, regulation of apoptosis). As a case study, we used the contextualized interactome to highlight the cellular processes induced upon Mycobacterium tuberculosis infection. CONCLUSION: Our work confirms that contextualizing interactomes improves the biological significance of bioinformatic analyses. More specifically, studying such inferred network rather than focusing at the gene expression level only, is informative on the processes involved in the host response. Indeed, important immune features such as apoptosis are solely highlighted when the spotlight is on the protein interaction level.